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^4) Procede de visualisation tridimensionnelle par lancer de rayons avec interpolation locale. 



(57) L'invention concerne tous les domaines qui requierent 
aes moyens de visualisation tridimensionnelle et dans les- 
quels les bases de donnees sont decrites a ('acquisition 
comme un empilement de coupes bidimensionnelles ou di- 
rectement en 3-D mais avec un echantillonnage anisotrope 
dans les trois dimensions. Le procede interesse done, en- 
tre autres, toutes les sources d'imagerie medicale fournis- 
sant ce type de donnees. 

Traditionnellement 1'isotropie des bases decrites ci- 
dessus est retablie, avant toute manipulation ou visualisa- 
tion, par interpolation de coupes. Ce procede presente Tin- 
convenient d'augmenter le volume des donnees et done 
d'accrottre la taille de la memo ire necessaire sur les syste- 
mes de traitement. 

L'invention se caracterise en ce qu'aucune interpolation 
globale prealable n'est realisee sur la base de donnees 
anisotrope (7) avant la visualisation de celle-ci par lancer 
de rayons. Pour eel a elle suit le processus suivant: 

- (8) lancer du rayon dans la base anisotrope, 

- (9) estimation de la position du point d'intersectton du 
rayon avec le sujet a visualiser dans la base de donnes, 

- M0) interpolation sur le voisinage de ce point estime, 

- (11) detection fine de la surface dans cette zone inter- 
polee, 

- (12) ombrage en fonction des parametres contenus 
dans cette zone. 
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Precede de visualisation tridimenslonnelle par lancer de rayons avec 
interpolation locale. 



L 1 invention concerne tous les domaines qui requierent des moyens 
de visualisation tridimenslonnelle (3-D) et dans lesquels les bases de 
donnees sont decrites a 1' acquisition comme un empilement de coupes 
bidimensionnelles ou directement en 3-D mais avec un echantillonnage 
10 anisotrope dans les trois dimensions. Le precede interesse done, entre 
autres, toutes les sources d* imagerie medicales fournissant ce type de 
donnees : II peut etre implante sur leur informatique hote mais aussi 
sur tout calculateur specialise ou d' usage general. 

15 Les techniques de visualisation tridimenslonnelle a ctuel lenient 

peuvent etre identifies par le type de representation des objets 
contenus dans la base de donnees qu' elles adoptent. II s'agit : 

- soit de representations surfaciques (description de la surface de 
l'objet par des approximations d'ordre 1 ou plus, ou encore par des 

20 voxels -volume elementaire generalement cubique) , 

- soit de representations volumiques, qui contrairement aux preceden- 
tes conservent V ensemble des informations presente dans la base de 
donnees initiale. 

25 Dans chacun de ces grands cadres, de nombreuses methodes ont ete 

proposees. L* invention concerne la representation par volume et plus 
particulierement les techniques par lancer de rayons (par opposition a 
des techniques de type projection ou de type codage par arbre octal ou 
"octree" ) . 
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Les approches basees sur le lancer de rayons sont apparues dans 
les annees 70 en synthese d' images : la recherche d'effets speciaux 
(reflexion multiple, refraction, ...) a cependant eu comme consequence 
de rendre complexes leurs applications. Leur utilisation pour le rendu 
de volumes ne remonte qu'a trois ou quatre ans . La visualisation de 
volumes apporte pourtant une simplification majeure, des lors que des 
perceptions realistes des bases de donnees seront privilegiees par 
rapport a la recherche d'effets speciaux. Seuls alors les rayons pri- 
maires seront a calculer. Ces techniques sont maintenant parfaitement 
eprouvees et dans le domaine medical on pourra se referer aux articles 
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de Drebin (DREBIN R. , CARPENTER L. et HANRAHAN P., Volume rendering, 
ACM Com. Graph., 12, 4, Aout 88, pp 65-74), de Levoy (LEVOY M. , 
Display of surfaces from volume data, Thesis, University of North 
Carolina at Chapel Hill, 1989), de Robb (ROBB R. and BARRILLOT C, 
Interactive 3-D display and analysis, SPIE vol 939 Hybrid Image and 
Signal Processing, 1988, pp 173-202) et Trousset (TROUSSET. Y. , 
Visualisation d'objets medicaux tridimensionnels , These, ENST Paris, 
sept 1987 et brevet francais FR 2.613.509). 

Tous ces travaux ont apporte des ameliorations au schema de base 
du lancer de rayons, soit sur le plan de la qualite des representa- 
tions, soit en terme de temps de calcul. La methode, decrite par la 
suite, a l'avantage de reduire les capacites memoires requises pour la 
representation 3-D et ceci sans limiter les autres f onctionnalites 
(transformations geometriques , decoupes par des plans quelconques , 
detection de sous obj ets , qualites de 1' image, ...). 

En visualisation d 1 image 3-D, un compromis doit etre trouve entre 
trois facteurs principaux : 

(a) la qualite des images et des fonctions offertes a 1 ' utilisateur , 

(b) les temps de reponse du systeme en sachant que 1 * interactivity 
homme-machine (meme pour des tailles elevees de bases de donnees) 
doit etre assuree, 

(c) les couts des stations d'imagerle. 

L' invention proposees permet d' apporter un gain significatif sur 
le trois ieme point (c) en diminuant notamment la taille de la memoirs 
necessaires aux systemes de traitement. 

La particularity de la plupart des bases de donnees 3-D 
constitutes a partir de coupes (ou plans) 2-D est de presenter une 
anisotropie axiale. C'est-a-dire que la resolution -1 ' echantillonnage- 
*dans le plan est meilleure que dans la direction orthogonale. 

Tous les travaux, publies k ce jours, retablissent avant toute 
visualisation ou manipulation de la base de donnees I'isotropie du 
volume en interpolant entre les coupes. La taille totale de la base de 
donnees est ainsi multipliee par un facteur N, avec N>2. 

II apparait avantageux, pour eviter cette augmentation du volume 
de donnees, de ne pas pratiquer une interpolation globale prealable et 
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de visualiser directement le contenu de la base anisotrope. Une seule 
publication (TROUSSET. Y. , Visualisation d'objets medicaux tridimen- 
sionnels, These, ENST Paris, sept 1987) propose un schema permettant 
de faire l'economie d'une telle augmentation mais en contraignant le 
probleme de telle sorte que la condition (a) n'est plus totalement 
respectee. Sa methode de visualisation d'une base de donnees tridimen- 
sionnelle formees de coupes empilees est caracterise par le fait que 
son incidence de vision est choisie de maniere a ce que I'ecran est 
place perpendiculairement aux coupes et que les rayons suivent une^ 
geometrie parallele. Ainsi, a chaque ligne de I'ecran correspond 
1' image de la tranche d'une coupe. Apres un lancer de rayons sur la 
base anisotrope, il lui suffit d f interpoler des ligne entre les lignes 
de I'ecran. Cette methode n'est generalisable ni dans le cas ou 
1* incidence de vision est quelconque, ni dans celui d'une geometrie 
conique . 

La presente invention vise a fournir un procede de visualisation 
par lancer de rayons ne necessitant pas d' interpolation prealable de 
la base de donnees et s ' af f ranchissant des restrictions citees 
prealablement . 

C'est-a-dire que : 

- la base de donnees anisotrope initiale reste inchangee, 

- seules les parties visibles (surfaces ou volumes internes) a 
l 1 observateur et independamment du point de we choisi devront etre 
interpolees lors de la visualisation. 

Dans ce but, 1' invention propose une methode caracterisee par une 
sequence d' operations a realiser pour chaque rayon et qui se decompose 
en differente phases : 

(1) definition du rayon (traj ectoire) dans la base anisotrope, 

(2) estimation de la position du point d* intersection du rayon avec le 
ou les sujets a visualiser dans la base de donnes , 

(3) interpolation sur le voisinage de ce point estime (le choix de 
1' interpolateur est libre, il faut seulement qu' il soit local), 

(4) recherche sur le rayon dans la zone interpolee du point 
d 1 intersection, 

(5) calcul eventuel de la normale, 

(6) specification de la valeur au plan image (ombrage) . 
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L' invention sera mieux comprise a la lecture de la description 
qui suit, d'un mode particulier de mise en oeuvre de l 1 invention, 
donne a titre d'exemple non limitatif. Cet exemple reprend les 
principes de base du lancer de rayons et definit les notations 
5 indispensables . Dans un but de simplification, l'exemple est decrit en 
geometrie parallele. Son principe s'etend cependant aux projections 
perspectives . 

La description se refere aux schemas qui 1 1 accompagnent dans 
lesquels : 

10 - la figure 1 est un schema qui presente les deux reperes (ecran 

et objets) et les relations geometriques entre eux, 

- la figure 2 est un schema qui represente le lancer de rayons en 
geometrie parallele, 

- la figure 3 est un graphe destine a illustrer la suite d'opera- 
15 tions du lancer de rayons avec interpolation locale, 

- la figure 4 est un schema qui pose le probleme de 1' intersec- 
tion entre le rayon et l'objet decrit dans la base de donne es , 

- la figure 5 est un schema qui illustre une des methodes d' esti- 
mation approchee de la position de 1' intersection rayon/objet, 

20 - la figure 6 est un schema qui illustre une deuxieme methode 

d* estimation approchee de cette position, 

- la figure 7 presente des histogrammes permettant 1' evaluation 
des methodes proposees. 

25 

Un rayon dans une base isotrope est definit de la maniere 
suivante : 

La scene 3-D possede son propre repere (1) (voir figure 1) appele 
repere objet R^. Celui-ci depend de la disposition des coupes 

30 tomographiques . Traditionnellement ces dernieres sont considerees 
paralleles au plan (x ob ,y ob ) et empilees selon l'axe z ob . 

L' autre systeme de coordonnees (2) -le repere ecran ou repere 
image R^- est lie a la visualisation. L 1 ecran represente le plan 
(x ec ,y ec ) . La scene est placee du cote des z ec positifs. 

35 Ces deux reperes sont lies par des relations de rotation et de 

translation. Ces relations dependent des pararaetres de vision. La 
direction de 1' observation est defini e par ses coordonnees spheriques 
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dans le repere objet (1) qui forment une matrice de rotation [R] entre 
les deux reperes . Les origines des reperes sont liees par un vecteur "? 
de translation. 

Lors du lancer d'un rayon, la position relative d'un point de ce 
dernier par rapport a 1* objet doit etre connue en permanence. Ceci 
impose constannnent de passer d'un repere a 1" autre. Toutefois, pour 
visualiser une base de donnees formee de voxels, il est plus efficace 
de tourner I'ecran autour de 1' objet que de tourner 1* objet lui-meme. 
Ce sont done les transformations ecran/objet qui interviendront le 
plus souvent . 

Ceci revient a appliquer la formule suivante : 

[Xob.Yob-Zob] = [R3- [x ec ,y ec ,z ec ]+T 

Le principe du lancer de rayons est le suivant (voir figure 2). A 
partir du point de depart (3) P 0 (x 0 ,y 0 , z 0 ) , le rayon (4) est parcouru 
de I'ecran (5) vers la scene (6). A chaque increment, une transforma- 
tion de repere ecran/objet est effectuee. Sur le resultat de celle-ci 
-c 1 est-a-dire les coordonnees de la position du rayon dans le repere 
objet- sont verifiees les conditions d' intersection (par une fonction 
que nous appellerons COND1TION_D_INTERSECTION) . Ces dernieres peuvent 
etre fondees sur la densite, le gradient local ou bien d'autres 
attributs du voxel. Une fois l 1 intersection constatee, la couleur du 
pixel (x 0 ,y 0 ) est calculee. 

Le rayon est une droite dans le repere ecran. Les parametres de 

ce rayon sont : 

- son point de depart Po (xq.vq.zq), 

- son vecteur directeur d [ 6^, dy, d z j , 

- une variable t representant la position le long du rayon. 
Son equation s ' ecrit : 



r x(t) - x 0 + d x .t 
<(y(t:) = y 0 + dy.t 
_z(t) - z 0 + d z .t 



Si les rayons sont paralleles entre eux et a 1* 
d^vaut [0,0,1]. 



axe z. 



■ec 
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L' equation du rayon se ramene a : fx(t) — xq 

< y(t) - y 0 

( z(t) - z 0 + t 

c'est a dire : [x ec (t) ,y ec (t) , z ec (t) ] = [x 0 ,y 0 ,z 0 +t] 
ou x 0 et y 0 sont les coordonnees du pixel de depart, . 

Les coordonnees du rayon dans le repere objet, pour un t donne , sont 
[x ob (t),y ob (t),z ob (t)] - [R].[x ec (t),y ec (t),z^ c (t)3 + T* 
- W ■ [x 0 ,y 0 ,zo +t ] + T 

Cette operation se decompose en trois termes : 



iXob(t) - x o-rii + y 0 .r 2 i + (z 0 +t).r 31 + t x 
y ob (t) - x 0 .r 12 + y 0 .r 22 + (z 0 +t).r 32 + t y 
ZobCt) = x 0- r 13 + yo- r 23 + (Z0 +t )- r 33 + c z 
avec r,-j les compos ante de la ma trice [R] et t| les valeurs du 

vecteur ?. 

La transformation des coordonnees & chaque pas du rayon ralenti 
20 le processus. Nous incrementerons done directement le rayon dans le 
repere objet. 

Ce L increment vaut : 

iXebCtf) " Wti-l) = (ti-ti-l) - r 31 
y 0 b(ti) - y 0 b(tf-i> = (ti-t i . 1 ).r 32 

Si nous considerons At = tf-ti.-j (pas du rayon), alors nous avons : 

+ [ r 31» r 32' r 33] -At 

30 



En rappel et en reference avec le schema synoptique de la figure 
3 voici les differentes operations du lancer de rayons sur une base 
35 anisotrope (7) (cette derniere, dans cet exemple , sera composee de 
coupes presentant une seule anisotropic selon z ob ) : 

(8) lancer du rayon dans la base anisotrope, 

(9) estimation de la position du point d' intersection du rayon avec le 
ou les sujets a visualiser dans la base de donnes , 
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(10) interpolation sur le voisinage de ce point estime, 

(11) detection fine de la surface dans cette zone interpolee, 

(12) ombrage en fonction des parametres contenus dans cette zone. 



Le rayon dans la base anisotrope est definit de la maniere 
suivante : 

L' echantillonnage du repere R^ n'est plus isotrope. Si D 
10 represente la distance inter-coupes et p la taille d'un pixel des 
coupes originales, le rapport de compression (note e) de R ob selon 
1 1 axe z 0 b vaut : e-D/p . 

Le re -echantillonnage isotrope, par interpolation de e-1 coupes 
f ictives entre deux coupes originales , definit un nouveau repere 
15 appele R fso . La transformation entre R iso et R ob est soumis a une 
matrice de mise a l'echelle E. 

10 0 
E - 0 1 0 
0 0 1/e 

20 Le lancer de rayons, decrit precedemment est effectue du repere 

R ec sur Rf S0 . Les transformations geometriques entre les differents 
reperes sont : 

tx iso ,yi so ,z fso ] = [R] - [x ec ,y ec ,z ec ]+? > 

[ x ob.y 0 b» z ob3 = [E] . [x is0) y is0) z iso ] 

25 

L' estimation de la position de l 1 intersection du rayon avec le ou 
les sujets a visualiser dans la base de donnees peut se derouler 
30 ainsi: 

Voir figure 4; si pour les plans originaux (13) les rayons (14) 
touchent bien l'objet (15), pour les plans " intermedia ire s " une surfa- 
ce encore inexistante doit etre detectee. 

1' interpolation sur un voisinage volumique est gourmande en temps 
35 de calcul. La reduction de ce voisinage impose une estimation de la 
position de la surface la plus precise possible. 

Le module de detection est le suivant : le point d* intersection 
estime est un voxel du repere R iso qui appartient a la trajectoire du 
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rayon partant de (x ec ,y ec ) et qui repond a un critere appele 
ESTIMATION. C'est-a-dire : 



P a ( xa iso »y a iso» z3 iso) — v3 f soC^iso » y iso » 



z a i S0 ) £ Rayon(x ec ,yec) 



5 



si ESTIMATION^ 1 



^iso.y a iso.z a iso)= v rai 



10 



P a (x,y,z) 

v iso( x fso »yfso > z iso) 
Rayon(x eCf y ec ) 



point de detection approchee de la surface 
voxel de R iso de position (x iso , y iso , z iso ) 
rayon lance du pixel (x ec ,y ec ) 



ESTIMATION (x,y , z) : module d' estimation de la position de la surface 

Plusieurs modeles de la fonction ESTIMATION sont etudies , qui 
font intervenir des criteres geometriques et les proprietes physiques 
des structures. 



- l er modele 

Ce modele, qui se refere a la figure 5, est fondee sur la discre- 
tisations originale de la base de donnees . Un lancer de rayons sur une 
base aniso trope equivaut, en fait, a un lancer de rayons (16) sur une 

20 base isotrope pour laquelle chaque coupe (17) retrouve son epaisseur 
originale (18). Tout ce passe comme si dans le repere R fso les voxels 
des plans intermediates prenaient la valeur des plans originaux. Nous 
pouvons done detecter les voxels (18) de la surface de maniere class i- 
que c'est-a-dire des qu'un rayon (16) les penetrent (19). 

25 La condition de detection approchee s'ecrit : 

ESTIMATION(x a fso ,y a isoi z a 1 - so )-vrai 

si CONDITIONJV INTERSECTION(x a ob ,y a ob , z a ob +0 , 5)=vrai 
avec [x a ob ,y a ob ,z a ob ] = [E] . [x a iso ,y a iso , z a Tso ) 

C0NDITI0N_D_INTERSECTI0N est la fonction citee precedemment qui 
30 verifie si le voxel appartient ou non a l'objet. Le voxel candidat a 
pour adresse x a ob , y a ob et la partie entiere de z a ob +0 , 5 . Cette 
translation de 0,5 permet de repartir 1* epaisseur de part et d' autre 
de la coupe originale. 

35 - 2^ modele 

Ce modele se refere a la figure 6. Le rayon (20), calcule en 
nombre reel, est entoure d'autres rayons formes par les coordonnees 



15 
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entieres du rayon initiale. Ces rayons, que nous pouvons appeler 
rayons auxiliaires, touchent obligatoirement les zones connues de 
l'objet. Les coordonnees de la position approchee de 1* intersection du 
rayon sont determinees par interpolation des positions d' intersections 
des rayons auxiliaires avec l'objet. 

Seule 1' interpolation selon z^ est determinante . Les rayons 
auxiliaires peuvent etre reduits aux nombre de deux : celui qui se 
trouve au-dessus (z entier superieur) et celui qui se trouve en- 
dessous (z entier inferieur) . 

L' estimation se deroule de la facon suivante : les intersections 
des rayons auxiliaires superieurs (21) et inferieurs (22) sont 
recherchees par CONDITIONJD ' INTERS ECTION(x ob ,y ob , z ob +l) pour le rayon 
auxiliaire superieur et par CONDITIONED' INTERS ECTI ON (x ob) y ob , z ob ) pour 
le rayon auxiliaire inferieur. Le point de la surface estime est 
represente par l 1 intersection (23) du rayon (20) et de la ligne (24) 
definie par les deux points d" intersection trouves (21) et (22). 

- 3^ modele 

Le premier modele peut etre affine. Le parametre physique 
(densite, ...) associe au point detecte est calcule par interpolation 
lineaire des valeurs des coupes environnantes . Si cette valeur repond 
aux conditions d* intersection, ce point est retenu, sinon le rayon est 
increments . 

ESTIMATION (x a iso , y a iso , z a iso ) =vrai 

si CONDITI0N_D , INTERSECTI0N(v 1 - so (x a 1 - so ,y a is01 z a 1 - s0 ))=vrai 
avec 

v iso( xa isot y a iso» za iso) = 

( ( e - d ) /e ) . v ob ( x a ob , y a ob , z a ob ) + ( d/e ) . v ob ( x a ob , y a ob , z a ob +l ) 

e : rapport de compression 

d : distance entre v a iso et v ob iso £ coupes originales selon z 
et [x%b>y a ob. za ob] = [EMx^ 0I y a iS0i z a f80 ] 
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Le module de 1 ' interpolation locale est le suivant : 
Toutes les techniques d' interpolation a partir de donnees locales 
sont utilisables qu'elles soient du premier ordre (interpolation 
lineaire) ou plus elaborees (fonctions splines, ... ). Cette liste 
n'est pas limitative, la seule restriction est le caractere local de 
l 1 interpolation. 

La taille et la forme de la zone environnant le point d' intersec- 
tion estime et dans laquelle sera menee 1 1 interpolation peut etre 
definie arbitrairement. Different criteres permettent toutefois un 
ajustement de ces parametres (cette liste est non limitative) : 

- la precision du modele de detection, 

- le rapport e (si e est faible -2 ou 3- la position estimee est 
proche de la surface "reelle"), 

- la technique d' interpolation retenue. II s'avere que les points 
estimes par l'operateur de detection sont proches de la surface de 
l'objet apres interpolation lineaire. Si la technique d' interpola- 
tion locale est lineaire, la region peut se resumer aux proches 
voisins (connexes) du voxel detecte. Pour les techniques d* ordre 
superieur cette region doit etre etendue , 

- la direction de vision. La region doit etre etendue selon la direc- 
tion de vision (celle du rayon) , 

y la methode de calcul de l'orabrage (des voisins qui y sont 
necessaires) . 

La detection fine de la surface et 1'ombrage se deroulent ainsi : 
Les differents voxels, issus du module d' interpolation, sont 
stockes dans une petite base auxiliaire. Le rayon, transpose dans 
cette base auxiliaire et isotrope, procede de maniere classique 
(recherche de la surface, ombrage -en fonction du gradient de la 
surface, de 1 ' eloignement de l'ecran, des sources de lumieres, de la 
nature de la surface, . . . - , . . . ) . 

La validation de la methode a ete effectuee sur deux plans : 

- 1* appreciation subjective des differences entre images par 
comparaison visuelle, 

- 1* analyse quantitative des valeurs issues des procedures proposees. 
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Cette derniere pouvait etre realisee directeraent dans l'espace 
3-D interpole apres detection. La valeur locale des comparaisons et la 
complexity de mise en oeuvre ont conduit a preferer une comparaison 
des images projetees sous differents points de vue. La procedure 
5 appliquee est la suivante : 1) acquisition de bases de donnees 3-D 
anisotrope; 2) interpolation globale et visualisation sous un angle 
d' incidence particulier; 3) pour le meme point de vue, calcul de 
1' image par interpolation locale sur le rayon; 4) soustraction point 
par point des images obtenues en (2) et (3) et statistiques sur les 
10 differences . 

Cette evaluation a ete faite pour des nombres de coupes 
interpolees variables, en utilisant les methodes de detection decrites 
precedemment et differents schemas d' interpolation. Les resultats 
apres reparcours local du rayon ne presentant pas de difference 
15 significatives, seuls les points de detection approchee ont ete 
utilises, ce qui a comme interet de simplifier la methode . 

Pour I'une des bases de donnees ayant fait l'objet de test, base 
de resolution 128 X 128 X 136 voxels (volume isotrope) , la suppression 
d'une coupe sur deux a ete effectuee. La base de donnees, rendue ainsi 
20 anisotrope, a ete visualisee soit apres interpolation globale 

(lineaire), soit par lancer de rayons avec interpolation locale. La 
figure 7 nous montre les histogrammes des differences entre 1' image 
originale (base de donnees isotrope de depart) et les images produites 
apres interpolation globale (25) ou interpolation locale (l er modele 
25 de detection (26), 2^ (27) et 3^ (28)). Elle montre que le 3^ 

modele de detection donne les meilleurs resultats et que ceux ci sont 
extremement proche de ceux de 1* interpolation globale. 

Naturellement, la presente invention n'est pas limitee a la 
30 description qui vient d'etre faite. En particulier elle peut s'appli- 
quer a toute les techniques d' acceleration du processus de lancer de 
rayons. L'avantage de ce procede est qu* il peut etre generalisable ou 
s'inserer a toute sequence d'operateurs le long du rayon exploitants 
les distribution locales des parametres physiques ou simplement des 
35 criteres geometriques . 
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Revendications 

1 - Procede de visualisation tridimensionnelle du contenu d'une 
base de donnees volumique presentant une ou plusieurs anisotropies 
axiales (echantillonnage different selon les trois axes) , selon une 
incidence quelconque, utilisant le principe de lancer de rayons et 
caracterise en ce qu 1 il comprend les operations suivantes : 

- definition du rayon (trajectoire) dans la base aniso trope, 

- estimation de la position du point d' intersection du rayon avec le 
ou les sujets k visualiser dans la base de donnees, 

- interpolation sur le voisinage de ce point estime (le choix de 
l 1 interpolation est libre, il faut seulement qu' il soit local), 

- recherche sur le rayon dans la zone interpolee du point 
d' intersection, 

- calcul eventuel de la normale en ce point et specification de la 
valeur au plan image (ombrage) . 

2 - Procede selon la revendication 1, caracterise en ce que la 
position du point d 1 intersection du rayon avec le ou les sujets 
contenus dans la base anisotrope est estimee lorsque le rayon penetre 
un voxel (anisotrope) de cette base, voxel qui repond aux conditions 
d' intersection souhaitees . 

3 - Procede selon la revendication 1, caracterise en ce que la 
position du point d' intersection du rayon avec le ou les sujets 
contenus dans la base anisotrope est estimee par une methode qui se 
sert du principe suivant : le rayon, calcule en nombre reel, est 
entoure d'autres rayons formes par les coordonnees entieres du rayon 
initiale, que ces rayons, que nous pouvons appeler rayons auxiliaires, 
touchent obligatoirement les zones connues de I'objet; ainsi les 
coordonnees de la position approchee de. 1 1 intersection du rayon sont 
determinees par interpolation des positions d 1 intersections des rayons 
auxiliaires avec I'objet. 
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4 - Procede selon la revendication 3, caracterise en ce que la 
position du point d' intersection du rayon avec le ou les sujets 
contenus dans la base anisotrope est estimee par une methode qui se 
sert du fait que si la base de donnees ne presente qu'une seule 

5 anisotropic axiale seule 1' interpolation selon cette direction est 
determinate; ainsi les rayons auxiliaires peuvent etre reduits au 
nombre de deux : ceux qui se trouvent de part et d* autre du rayon 
selon cette direction; 1' estimation se deroule le principe suivant ; 
les intersections des deux rayons auxiliaires sont recherchees , le 

10 point de la surface estime est represents par 1* intersection (23) du 
rayon (20) et de la ligne (24) definie par les deux points 
d 1 intersection trouves (21) et (22). 

5 - Procede selon la revendication 2, selon lequel le parametre 
15 physique associe au point trouve selon le procede de la revendication 

2 est calcule par interpolation lineaire des valeurs des coupes 
environnantes ; la position du point d' intersection du rayon avec le ou 
les sujets contenus dans la base anisotrope est estimee si la valeur 
du parametre physique repond aux conditions d' intersection. 



20 



25 



6 - Procede selon i'une quelconque des revendications 1 a 5, 
caracterise en ce que la position estimee de 1 ' intersection entre le 
rayon et le ou les sujets contenus dans la base de donnees est 
considere comme position exacte de cette intersection. 



7 - Procede selon I'une quelconque des revendications la 6, 
caracterise en ce que la taille et la forme de la zone environnant le 
point d» intersection estime et dans laquelle sera menee 1' interpola- 
tion sont definies soit arbitrairement soit selon un ou plusieurs 
criteres de la liste suivante : 
30 - la precision du modele de detection, 

- le ou les rapports caracterisant 1' anisotropic , 

- la technique d* interpolation retenue, 

- la direction de vision, 

- la methode de calcul de I'ombrage (des voisins qui y sont 
35 necessaires) . 
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8 - Procede selon l'une quelconque des revendications 1 a 7, 
caracterise en ce que les differents voxels, issus du module 
d' interpolation, sont stockes dans une petite base auxiliaire 
temporaire et que le rayon, transpose dans cette base auxiliaire et 
isotrope, procede de maniere classique pour trouver le ou les objets a 
visualiser et pour specifier la valeur du plan image (ombrage) . 
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Figure 1 
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PLRNCHE 3/7 
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Definition du rayon 
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Interpolation sur 1e voisinage 
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Detection fine de U surface 
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Figure 3 




BNSDOCID: <FR_2MSS2AMJ_?- 




BNSDOaD: <FR__2e82S24A1JU» 




BN8DOCID: <PR__2ee2624A1JU» 



REPUBLIQUE FRAN^aTSE 



INS TIT UT NATIONAL 
de la 

PROPRIETE INDUSTRIELLE 



2662524 



N° d'enregistrement 
national 



RAPPORT DE RECHERCHE 

etabli sur la base des dernier es revend teat ions 
deposees avant le commencement de la recherche 



FR 9006683 
FA 443938 



Categorie 



DOCUMENTS CONSIDERES COMME PERTINENTS 



Citation du document avec indication, en cas de besoin, 
des parties pertinentes 



Revindications 
concernees 
de la demande 
examinee 



SYSTEMS & COMPUTERS IN JAPAN, vol. 18, 
no. 11, novembre 1987, pages 31-39, 
Silver Spring, MD, US; TAKAAKI AKIMOTO 
et al . : "Pixel-selected ray-tracing" 
* Page 41, colonne 1, lignes 7-10; page 

41, colonne 2, lignes 10-13,21-24; page 

42, colonne 1, lignes 

1-4,19-22,26-28,35-36; page 43, colonne 
1, lignes 23-28; page 43, colonne 2, 
lignes 51-54; page 44, colonne 1, 
lignes 1-2,18-19; page 44, colonne 1, 
ligne 26; page 45, colonne 1, lignes 
18-24; page 47, colonne 2, lignes 19-24 



EP-A-0 362 123 (B.T. LEE) 
* Colonne 1, lignes 29-40; colonne 1, 
ligne 47 - colonne 2, ligne 2; colonne 
2, lignes 40-52 * 



1-7 



1-2,7-8 



DOMAEVES TECHNIQUES 
RECHERCHES (\au CL5) 



G 06 F ■ 15/72 



Date <T zebivement de ta recherche 

05-02-1991 



Exanrinatenr 

BARBA M. 



CATEGORIE DES DOCUMENTS CITES 

X : particulierement pertinent i lui seul 

Y : particulierement pertinent en combinaison avec un 

autre document de la meme categorie 
A : pertinent a 1'encontre d'au moms une revendication 

au arriere-plan technologique general 
O : divulgation non-ecrite 

P : document intercalaire 



T : theorie ou principe i la base de 1'fnventlon 

E : document de brevet beneficiant d'ime date anterieure 

a la date de depot et qui n'a etc publfe qu'a cette date 

de depfit ou qu'a une date posterieure. 
D : cite dans la demande 
L : cite pour d'autres ralsons 

& : raembre de la meme famille, document correspondant 



BNSOOC1D: <FR__29«2S24A1JL> 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 



□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




BLURRED OR ILLEGIBLE TEXT OR DRAWING 



^ADED TEXT OR DRAWING 



